High-throughput sequencing of insect specimens with sub-optimal DNA preservation using a practical, plate-based Illumina-compatible Tn5 transposase library preparation method.

Entomological sampling and storage conditions often prioritise efficiency, practicality and conservation of morphological characteristics, and may therefore be suboptimal for DNA preservation. This practice can impact downstream molecular applications, such as the generation of high-throughput genomic libraries, which often requires substantial DNA input amounts. Here, we use a practical Tn5 transposase tagmentation-based library preparation method optimised for 96-well plates and low yield DNA extracts from insect legs that were stored under sub-optimal conditions for DNA preservation. The samples were kept in field vehicles for extended periods of time, before long-term storage in ethanol in the freezer, or dry at room temperature. By reducing DNA input to 6ng, more samples with sub-optimal DNA yields could be processed. We matched this low DNA input with a 6-fold dilution of a commercially available tagmentation enzyme, significantly reducing library preparation costs. Costs and workload were further suppressed by direct post-amplification pooling of individual libraries. We generated medium coverage (>3-fold) genomes for 88 out of 90 specimens, with an average of approximately 10-fold coverage. While samples stored in ethanol yielded significantly less DNA compared to those which were stored dry, these samples had superior sequencing statistics, with longer sequencing reads and higher rates of endogenous DNA. Furthermore, we find that the efficiency of tagmentation-based library preparation can be improved by a thorough post-amplification bead clean-up which selects against both short and large DNA fragments. By opening opportunities for the use of sub-optimally preserved, low yield DNA extracts, we broaden the scope of whole genome studies of insect specimens. We therefore expect these results and this protocol to be valuable for a range of applications in the field of entomology.

Platelet and mitochondrial RNA is decreased in plasma-derived extracellular vesicles in women with preeclampsia-an exploratory study.

BACKGROUND: Circulating extracellular vesicles (EVs) are increased in preeclampsia (PE) and are associated with severity and progression. We examined in this exploratory cohort study if the mRNAs and long noncoding RNAs (lncRNAs) in plasma-derived EVs were dysregulated in PE compared to normal pregnancy and display different temporal patterns during gestation. METHODS: We isolated EVs from plasma at weeks 22-24 and 36-38 in women with and without PE (n=7 in each group) and performed RNA-seq, focusing on mRNAs and lncRNAs. We validated highly expressed mitochondrial and platelet-derived RNAs discovered from central pathways in 60 women with/without PE. We examined further one of the regulated RNAs, noncoding mitochondrially encoded tRNA alanine (MT-TA), in leukocytes and plasma to investigate its biomarker potential and association with clinical markers of PE. RESULTS: We found abundant levels of platelet-derived and mitochondrial RNAs in EVs. Expression of these RNAs were decreased and lncRNAs increased in EVs from PE compared to without PE. These findings were further validated by qPCR for mitochondrial RNAs MT-TA, MT-ND2, MT-CYB and platelet-derived RNAs PPBP, PF4, CLU in EVs. Decreased expression of mitochondrial tRNA MT-TA in leukocytes at 22-24 weeks was strongly associated with the subsequent development of PE. CONCLUSIONS: Platelet-derived and mitochondrial RNA were highly expressed in plasma EVs and were decreased in EVs isolated from women with PE compared to without PE. LncRNAs were mostly increased in PE. The MT-TA in leukocytes may be a useful biomarker for prediction and/or early detection of PE.

A 39 kb structural variant causing Lynch Syndrome detected by optical genome mapping and nanopore sequencing.

Lynch Syndrome (LS) is a hereditary cancer syndrome caused by pathogenic germline variants in one of the four mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2. It is characterized by a significantly increased risk of multiple cancer types, particularly colorectal and endometrial cancer, with autosomal dominant inheritance. Access to precise and sensitive methods for genetic testing is important, as early detection and prevention of cancer is possible when the variant is known. We present here two unrelated Norwegian families with family histories strongly suggestive of LS, where immunohistochemical and microsatellite instability analyses indicated presence of a pathogenic variant in MSH2, but targeted exon sequencing and multiplex ligation-dependent probe amplification (MLPA) were negative. Using Bionano optical genome mapping, we detected a 39 kb insertion in the MSH2 gene. Precise mapping of the insertion breakpoints and inserted sequence was performed by low-coverage whole-genome sequencing with an Oxford Nanopore MinION. The same variant was present in both families, and later found in other families from the same region of Norway, indicative of a founder event. To our knowledge, this is the first diagnosis of LS caused by a structural variant using these technologies. We suggest that structural variant detection be performed when LS is suspected but not confirmed with first-tier standard genetic testing.

Vitamin D in Alzheimer's Disease: Low Levels in Cerebrospinal Fluid Despite Normal Amounts in Serum.

BACKGROUND: Vitamin D insufficiency has been suggested as a dementia risk factor. OBJECTIVE: In this cross-sectional, explorative study we investigated whether levels of vitamin D in cerebrospinal fluid (CSF) are lower in patients with positive biomarkers of Alzheimer's disease (AD) compared to cognitively healthy controls and whether polymorphisms of the vitamin D receptor (VDR) gene, FokI, BsmI, ApaI, and TaqI, are associated with levels of vitamin D in CSF and cognition. METHODS: We included 100 patients≥65 years assessed for cognitive impairment and 76 cognitively healthy controls. Levels of 25-hydroxyvitamin D (25(OH)D) in both serum and CSF, and VDR polymorphisms were analyzed. RESULTS: The mean level of 25(OH)D in serum was 78.6 (SD 28.9) nmol/l. While serum levels of 25(OH)D were not significantly different between the groups, CSF levels of 25(OH)D were significantly lower in patients with positive AD core biomarkers (p = 0.001) compared to patients without such biomarkers. Individuals with the BsmI major homozygote genotype had significantly lower results on a 10-word delayed recall test (p = 0.044) and verbal fluency test (p = 0.013), and individuals with the TaqI major homozygote genotype had significantly lower results on a verbal fluency test (p = 0.030) compared to individuals with the corresponding minor homozygote genotype. CONCLUSION: Patients with positive AD core biomarkers have low CSF levels of 25(OH)D, despite sufficient serum levels. CSF levels of 25(OH)D do not seem to be affected by any of the four VDR gene polymorphisms. TaqI and BsmI major homozygote genotypes might be at increased risk for development of cognitive decline.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing.

BACKGROUND: Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 µg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits. RESULTS: We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits' fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection. CONCLUSIONS: The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Mapping the expression of an ANK3 isoform associated with bipolar disorder in the human brain.

The gene ankyrin-3 (ANK3) has been consistently associated with bipolar disorder (BD) in several genome-wide association studies (GWAS). The exact molecular mechanisms underlying this genetic association remain unknown. The discovery of a loss-of-function variant (rs41283526*G) in an alternatively spliced exon (ENSE00001786716) with a protective effect, suggested that elevated expression of this particular isoform could be a risk factor for developing the disorder. We developed a novel approach for measuring the expression level of all splice forms at a challenging genetic locus using a combination of droplet digital PCR and high-throughput sequencing of indexed PCR amplicons. The combined method was performed on a large collection of 568 postmortem brain samples of BD and SCZ cases and controls. We also studied the expression of the splice forms in a child-development cohort of 41 healthy males. We found that our approach can quantify the splice forms in brain samples, although with less precision than ddPCR. We detected highly significant differences in expression of splice forms and transcription start sites between brain regions, notably with higher expression of the BD-associated isoform in the corpus callosum compared to frontal tissue (mean fold change = 1.80, p < 1e-4). Although the patients in our sample expressed the BD-associated splice form at a similar level to controls, adolescents in our child-development cohort had a clearly higher expression level than younger children (mean fold change = 1.97, p = 5e-3). These results suggest that this ANK3 splice form may play a role in the myelin maturation of the human brain.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns.

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Rapid SARS-CoV-2 variant monitoring using PCR confirmed by whole genome sequencing in a high-volume diagnostic laboratory.

OBJECTIVES: The emerging SARS-CoV-2 variants of concern (VoC), B.1.1.7, B.1.351 and P.1, with increased transmission and/or immune evasion, emphasise the need for broad and rapid variant monitoring. Our high-volume laboratory introduced a PCR variant assay (Variant PCR) in January 2021 based on the protocol by Vogels et al. STUDY DESIGN: To assess whether Variant PCR could be used for rapid B.1.1.7, B.1.351 and P.1 screening, all positive SARS-CoV-2 airway samples were prospectively tested in parallel using both the Variant PCR and whole genome sequencing (WGS). RESULTS: In total 1,642 SARS-CoV-2 positive samples from individual patients were tested within a time span of 4 weeks. For all samples with valid results from both Variant PCR and WGS, no VoC was missed by Variant PCR (totalling 399 VoC detected). Conversely, all of the samples identified as "other lineages" (i.e., "non-VoC lineages") by the Variant PCR, were confirmed by WGS. CONCLUSIONS: The Variant PCR based on the protocol by Vogels et al., is an effective method for rapid screening for VoC, applicable for most diagnostic laboratories within a pandemic setting. WGS is still required to confirm the identity of certain variants and for continuous surveillance of emerging VoC.

An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Going low to reach high: Small-scale ChIP-seq maps new terrain.

Chromatin immunoprecipitation (ChIP) enables mapping of specific histone modifications or chromatin-associated factors in the genome and represents a powerful tool in the study of chromatin and genome regulation. Importantly, recent technological advances that couple ChIP with whole-genome high-throughput sequencing (ChIP-seq) now allow the mapping of chromatin factors throughout the genome. However, the requirement for large amounts of ChIP-seq input material has long made it challenging to assess chromatin profiles of cell types only available in limited numbers. For many cell types, it is not feasible to reach high numbers when collecting them as homogeneous cell populations in vivo. Nonetheless, it is an advantage to work with pure cell populations to reach robust biological conclusions. Here, we review (a) how ChIP protocols have been scaled down for use with as little as a few hundred cells; (b) which considerations to be aware of when preparing small-scale ChIP-seq and analyzing data; and (c) the potential of small-scale ChIP-seq datasets for elucidating chromatin dynamics in various biological systems, including some examples such as oocyte maturation and preimplantation embryo development. This article is categorized under: Laboratory Methods and Technologies > Genetic/Genomic Methods Developmental Biology > Developmental Processes in Health and Disease Biological Mechanisms > Cell Fates.

Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

Methylation Warfare: Interaction of Pneumococcal Bacteriophages with Their Host.

Virus-host interactions are regulated by complex coevolutionary dynamics. In Streptococcus pneumoniae, phase-variable type I restriction-modification (R-M) systems are part of the core genome. We hypothesized that the ability of the R-M systems to switch between six target DNA specificities also has a key role in preventing the spread of bacteriophages. Using the streptococcal temperate bacteriophage SpSL1, we show that the variants of both the SpnIII and SpnIV R-M systems are able to restrict invading bacteriophage with an efficiency approximately proportional to the number of target sites in the bacteriophage genome. In addition to restriction of lytic replication, SpnIII also led to abortive infection in the majority of host cells. During lytic infection, transcriptional analysis found evidence of phage-host interaction through the strong upregulation of the nrdR nucleotide biosynthesis regulon. During lysogeny, the phage had less of an effect on host gene regulation. This research demonstrates a novel combined bacteriophage restriction and abortive infection mechanism, highlighting the importance that the phase-variable type I R-M systems have in the multifunctional defense against bacteriophage infection in the respiratory pathogen S. pneumoniaeIMPORTANCE With antimicrobial drug resistance becoming an increasing burden on human health, much attention has been focused on the potential use of bacteriophages and their enzymes as therapeutics. However, the investigations into the physiology of the complex interactions of bacteriophages with their hosts have attracted far less attention, in comparison. This work describes the molecular characterization of the infectious cycle of a bacteriophage in the important human pathogen Streptococcus pneumoniae and explores the intricate relationship between phase-variable host defense mechanisms and the virus. This is the first report showing how a phase-variable type I restriction-modification system is involved in bacteriophage restriction while it also provides an additional level of infection control through abortive infection.

Pre-marked chromatin and transcription factor co-binding shape the pioneering activity of Foxa2.

Pioneer transcription factors (PTF) can recognize their binding sites on nucleosomal DNA and trigger chromatin opening for recruitment of other non-pioneer transcription factors. However, critical properties of PTFs are still poorly understood, such as how these transcription factors selectively recognize cell type-specific binding sites and under which conditions they can initiate chromatin remodelling. Here we show that early endoderm binding sites of the paradigm PTF Foxa2 are epigenetically primed by low levels of active chromatin modifications in embryonic stem cells (ESC). Priming of these binding sites is supported by preferential recruitment of Foxa2 to endoderm binding sites compared to lineage-inappropriate binding sites, when ectopically expressed in ESCs. We further show that binding of Foxa2 is required for chromatin opening during endoderm differentiation. However, increased chromatin accessibility was only detected on binding sites which are synergistically bound with other endoderm transcription factors. Thus, our data suggest that binding site selection of PTFs is directed by the chromatin environment and that chromatin opening requires collaboration of PTFs with additional transcription factors.

Monocytes accumulate in the airways of children with fatal asthma.

BACKGROUND: Activated T helper type 2 (Th2) cells are believed to play a pivotal role in allergic airway inflammation, but which cells attract and activate Th2 cells locally have not been fully determined. Recently, it was shown in an experimental human model of allergic rhinitis (AR) that activated monocytes rapidly accumulate in the nasal mucosa after local allergen challenge, where they promote recruitment of Th2 cells and eosinophils. OBJECTIVE: To investigate whether monocytes are recruited to the lungs in paediatric asthma. METHODS: Tissue samples obtained from children and adolescents with fatal asthma attack (n = 12), age-matched non-atopic controls (n = 9) and allergen-challenged AR patients (n = 8) were subjected to in situ immunostaining. RESULTS: Monocytes, identified as CD68+S100A8/A9+ cells, were significantly increased in the lower airway mucosa and in the alveoli of fatal asthma patients compared with control individuals. Interestingly, cellular aggregates containing CD68+S100A8/A9+ monocytes obstructing the lumen of bronchioles were found in asthmatics (8 out of 12) but not in controls. Analysing tissue specimens from challenged AR patients, we confirmed that co-staining with CD68 and S100A8/A9 was a valid method to identify recently recruited monocytes. We also showed that the vast majority of accumulating monocytes both in the lungs and in the nasal mucosa expressed matrix metalloproteinase 10, suggesting that this protein may be involved in their migration within the tissue. CONCLUSIONS AND CLINICAL RELEVANCE: Monocytes accumulated in the lungs of children and adolescents with fatal asthma attack. This finding strongly suggests that monocytes are directly involved in the immunopathology of asthma and that these pro-inflammatory cells are potential targets for therapy.

Transcriptomics of the Vaccine Immune Response: Priming With Adjuvant Modulates Recall Innate Responses After Boosting.

Transcriptomic profiling of the immune response induced by vaccine adjuvants is of critical importance for the rational design of vaccination strategies. In this study, transcriptomics was employed to profile the effect of the vaccine adjuvant used for priming on the immune response following re-exposure to the vaccine antigen alone. Mice were primed with the chimeric vaccine antigen H56 of Mycobacterium tuberculosis administered alone or with the CAF01 adjuvant and boosted with the antigen alone. mRNA sequencing was performed on blood samples collected 1, 2, and 7 days after priming and after boosting. Gene expression analysis at day 2 after priming showed that the CAF01 adjuvanted vaccine induced a stronger upregulation of the innate immunity modules compared with the unadjuvanted formulation. The immunostimulant effect of the CAF01 adjuvant, used in the primary immunization, was clearly seen after a booster immunization with a low dose of antigen alone. One day after boost, we observed a strong upregulation of multiple genes in blood of mice primed with H56 + CAF01 compared with mice primed with the H56 alone. In particular, blood transcription modules related to innate immune response, such as monocyte and neutrophil recruitment, activation of antigen-presenting cells, and interferon response were activated. Seven days after boost, differential expression of innate response genes faded while a moderate differential expression of T cell activation modules was appreciable. Indeed, immunological analysis showed a higher frequency of H56-specific CD4+ T cells and germinal center B cells in draining lymph nodes, a strong H56-specific humoral response and a higher frequency of antibody-secreting cells in spleen of mice primed with H56 + CAF01. Taken together, these data indicate that the adjuvant used for priming strongly reprograms the immune response that, upon boosting, results in a stronger recall innate response essential for shaping the downstream adaptive response.

Long Noncoding RNAs in Yeast Cells and Differentiated Subpopulations of Yeast Colonies and Biofilms.

We summarize current knowledge regarding regulatory functions of long noncoding RNAs (lncRNAs) in yeast, with emphasis on lncRNAs identified recently in yeast colonies and biofilms. Potential regulatory functions of these lncRNAs in differentiated cells of domesticated colonies adapted to plentiful conditions versus yeast colony biofilms are discussed. We show that specific cell types differ in their complements of lncRNA, that this complement changes over time in differentiating upper cells, and that these lncRNAs target diverse functional categories of genes in different cell subpopulations and specific colony types.

Transcriptome Remodeling of Differentiated Cells during Chronological Ageing of Yeast Colonies: New Insights into Metabolic Differentiation.

We present the spatiotemporal metabolic differentiation of yeast cell subpopulations from upper, lower, and margin regions of colonies of different ages, based on comprehensive transcriptomic analysis. Furthermore, the analysis was extended to include smaller cell subpopulations identified previously by microscopy within fully differentiated U and L cells of aged colonies. New data from RNA-seq provides both spatial and temporal information on cell metabolic reprogramming during colony ageing and shows that cells at marginal positions are similar to upper cells, but both these cell types are metabolically distinct from cells localized to lower colony regions. As colonies age, dramatic metabolic reprogramming occurs in cells of upper regions, while changes in margin and lower cells are less prominent. Interestingly, whereas clear expression differences were identified between two L cell subpopulations, U cells (which adopt metabolic profiles, similar to those of tumor cells) form a more homogeneous cell population. The data identified crucial metabolic reprogramming events that arise de novo during colony ageing and are linked to U and L cell colony differentiation and support a role for mitochondria in this differentiation process.

An optimized FAIRE procedure for low cell numbers in yeast.

We report an optimized low-input FAIRE-seq (Formaldehyde-Assisted Isolation of Regulatory Elements-sequencing) procedure to assay chromatin accessibility from limited amounts of yeast cells. We demonstrate that the method performs well on as little as 4 mg of cells scraped directly from a few colonies. Sensitivity, specificity and reproducibility of the scaled-down method are comparable with those of regular, higher input amounts, and allow the use of 100-fold fewer cells than existing procedures. The method enables epigenetic analysis of chromatin structure without the need for cell multiplication of exponentially growing cells in liquid culture, thus opening the possibility of studying colony cell subpopulations, or those that can be isolated directly from environmental samples.

Microbial Community Composition of Tap Water and Biofilms Treated with or without Copper-Silver Ionization.

Copper-silver ionization (CSI) is an in-house water disinfection method primarily installed to eradicate Legionella bacteria from drinking water distribution systems (DWDS). Its effect on the abundance of culturable Legionella and Legionella infections has been documented in several studies. However, the effect of CSI on other bacteria in DWDS is largely unknown. To investigate these effects, we characterized drinking water and biofilm communities in a hospital using CSI, in a neighboring building without CSI, and in treated drinking water at the local water treatment plant. We used 16S rDNA amplicon sequencing and Legionella culturing. The sequencing results revealed three distinct water groups: (1) cold-water samples (no CSI), (2) warm-water samples at the research institute (no CSI), and (3) warm-water samples at the hospital (after CSI; ANOSIM, p < 0.001). Differences between the biofilm communities exposed and not exposed to CSI were less clear (ANOSIM, p = 0.022). No Legionella were cultured, but limited numbers of Legionella sequences were recovered from all 25 water samples (0.2-1.4% relative abundance). The clustering pattern indicated local selection of Legionella types (Kruskal-Wallis, p < 0.001). Furthermore, one unclassified Betaproteobacteria OTU was highly enriched in CSI-treated warm water samples at the hospital (Kruskal-Wallis, p < 0.001).